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Abstract. Detailed simultaneous radiative transfer-statistical 
equilibrium calculations have been performed for lithium in 
atmospheres of solar metallicity G-M dwarfs and subgiants. 
Model atmospheres of stars with 3500 < T cB < 6000 K, 
3.0 < \ogg < 4.5 have been considered. The behaviour of 
level populations, source functions and curves of growth for 
the lithium doublet at A670.8 nm and two subordinate lines 
of astrophysical interest (A610.3 nm and A812.6 nm) are dis- 
cussed. For the coolest atmospheres of our grid we find that 
the dependence of non-LTE equivalent widths and profiles on 
gravity is very weak. The computed LTE and non-LTE curves 
of growth are given in tabular form. 
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1. Introduction 

The importance of lithium in many fields of astrophysics has 
been widely recognized during the last decades. The proceed- 
ings of a recent ESO/EIPC workshop (Crane 1995) describe 
the status of the present knowledge on this subject. The abun- 
dance of lithium is usually determined from the resonance dou- 
blet at A670.8 nm. The equivalent width of this feature, quite 
small in general, can reach values up to hundreds of mA in 
spectra of young pre-main sequence objects, where the dou- 
blet can get easily saturated. We restricted our study to stars 
with T c ff > 3500 K. Most lithium in the atmosphere of these 
stars is in the form of Li II ions. Therefore, the amount of neu- 
tral lithium can be sensitive to non-LTE (non local thermody- 
namic equilibrium) effects, notably overionization, leading to 
non-LTE effects on abundance determinations, even for weak 
lines. Strong lines in addition will have non-LTE affecting the 
source function. Studies of the non-LTE formation of Li I lines 
are required in order to get accurate abundances from the ob- 
served equivalent widths. 

We have presented our results on the non-LTE formation 
of lithium lines in the atmosphere of cool dwarfs and sub- 
giants in several papers (Magazzu et al. 1992, Martin et al. 
1994, Pavlenko 1995, Pavlenko et al. 1995), aimed mainly to 
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determine the lithium abundance in the atmosphere of pre- 
main sequence objects. Recently, Carlsson et al. (1994) pub- 
lished results of extensive non-LTE computations for lithium 
for a grid of model atmospheres with T e ff between 4500 and 
7500 K. They provide numerical corrections to get non-LTE 
Li abundances from LTE (local thermodynamic equilibrium) 
curves of growth. 

In this work we present new results of non-LTE studies 
concerning the Li I resonance line at A670.8 nm and also the 
subordinate lines at A610.3 nm and A812.6 nm, whose astro- 
physical interest has been recently recognized (Duncan 1991, 
Magazzu et al. 1995). The curves of growth presented here for 
the A670.8 nm line are a refinement of those used by Martin 
et al. (1994). 



2. The Procedure 

To solve the system of statistical balance equations and ra- 
diative transfer equations (the non-LTE problem) we followed 
the modification of the linearization method proposed by Auer 
and Heasley (1976). Our computation procedure of the non- 
LTE problem solution was described extensively in Pavlenko 
et al. (1995). Some details are given below: 

— we used a 20-level model of the Li I atom; 

— the 70 radiative and all possible collisional transitions were 
included into the rate matrix; 

— the ionization-dissociation equilibrium was computed for 
98 atoms and molecules including 6 lithium-containing 
molecules; 

— the multiplet radiative transitions were replaced by single 
transitions, even though absorption coefficient profiles were 
computed taking into account the multiplet structure; 

— the opacity due to atomic lines from Kurucz (1979) was 
taken into account. Additionally, we included the absorp- 
tion due to molecular bands in JOLA approximation, 
which gives good results for saturated bands (Tsuji 1994). 
Twenty-one molecular bands have been taken into account. 
For unsaturated molecular bands (in the case of stars with 
Tcft > 4000 K) results are not so good, but in this case 
the visible and violet parts of the spectrum are dominated 
by absorption due to atomic lines. All these opacities were 
treated in LTE; 
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— we used the correction factor E = 2 for the Unsold approx- 
imation of van der Waals damping for lithium lines (see 
Section 2.5 in Pavlenko et al. 1995); 

— we excluded from our computations the inelastic collisions 
of Li I atoms with hydrogen neutral atoms, due to the 
strong criticism raised (see Lambert 1993, Carlsson et al. 
1994) on Drawin's (1968, 1969) formulae and their modifi- 
cations by Steenbock and Holveger (1984). 

The main difference between this work and computations 
by Pavlenko et al. (1995) is the model atmospheres. Here we 
use models calculated with ATLAS9 and taken from Kurucz 
(1993) CD-ROM 13. We consider only solar metallicity mod- 
els, with microturbulence of 2 km s _1 . Any contribution from 
6 Li to the equivalent widths of the lines in study has been 
neglected. 

3. Results 

We treat here the range 3500 < T cff < 6000 K, which in- 
cludes atmospheres where the molecular absorption is relevant. 
In such atmospheres the thermalization processes increase due 
to the strong absorption by atomic and molecular lines in the 
frequencies of Li I resonance and subordinate lines. These pro- 
cesses are currently studied and will be discussed in forth- 
coming papers. Here we give a few results, relevant for un- 
derstanding the nature of non-LTE effects in Li I lines in cool 
atmospheres. For this purpose it is useful to consider the de- 
pendence of departure coefficients and source functions of the 
absorption lines on the optical depth. 
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3.1. The departure coefficients 

In Fig. [l| we present the departure coefficients of lithium levels 
bi = rii/n* versus ri M , the continuum optical depth at 1 /jm, 
in the atmosphere of a dwarf 4000/3.0, i.e. with T e g = 4000 K, 
logg = 3.0. Here n* and rn axe the LTE and non-LTE pop- 
ulations of the lithium levels, respectively. The model atom 
used was described in detail by Pavlenlso (1994). We present 
results of computations for abundanceail log iV(Li) = 0.5 and 
log 7V(Li) = 3.5, i.e. in the case of weak and strong lithium 
resonance doublet, respectively. 
We find that: 

— the behaviour of the departure coefficients depends on the 
strength of Li I resonance doublet, i.e. on the lithium abun- 
dance (see also Magazzu et al. 1992); 

— properly speaking, we should use the definition "overion- 
ization of lithium" to describe only the relevant process of 
formation of the statistical equilibrium of lithium. In fact, 
the statistical equilibrium of lithium depends on a whole 
set of radiative and collisional transitions (see Carlsson 
et al. 1994). Still, in the formation region of weak lithium 
A670.8nm lines (r = Tip ~ 10~ 2 ) overionization domi- 
nates. As a result we have here 6i < &2 < 1.0 (see also 
Steenbock & Holweger 1984); 

— in the case of saturated lithium resonance lines we have bal- 
ance of the radiative transitions (2s- 2p) in a large part of 
the stellar atmosphere. The population of the first lithium 
level is formed by chains of transitions from upper levels 
and continuum. In the formation region of strong resonance 

Throughout this paper log iV(Li) is given in the scale 
logiV(H) = 12 



Fig. 1. Departure coefficients of Li I levels in atmosphere of a 
4000/3.0 subgiant with a lithium abundance logA r (Li) = 0.5 (a) 
and logJV(Li) = 3.5 (b) 

lithium lines (r ~ 10~ J ) we have &i > 62. In the outer part 
of stellar atmospheres the first level may be even overpopu- 
lated with respect to LTE (see also Pavlenko 1992, Carlsson 
et al. 1994); 

— in addition, the departure coefficients of upper levels reach 
values greater than 1 in some regions of the atmosphere. 
These levels are more tied with continuum than with lower 
levels. 

3.2. The source function 

The nature and intensity of the non-LTE effects depend on 
Si/B v , i.e. the ratio of the source function to the Planck func- 
tion at the formation depth of a line. Here we study such a 
dependence for the lines considered in this work. In Fig. |^ we 
present the ratio Si/B„ computed for the Li I A670.8 nm res- 
onance doublet plotted against ri M for several model atmo- 
spheres with different effective temperatures and gravities. 
We find that in the A670.8 nm line formation region: 

— the ratio Si/B v is lower for stars with lower T e g; 

— the ratio is lower in atmospheres of lower luminosity stars; 

— the differences in Si/B v computed for logg = 3.0 and 4.5 
are very small for T e « = 3500 K. 

When the effective temperature decreases the lithium reso- 
nance lines (for a given lithium abundance) become stronger, so 
their thermalization depths shift towards the outer boundary 
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Fig. 2. Ratio Si/B v for the Li I resonance doublet in atmospheres 
5000/4.5 and 5000/3.0 (1), 4000/4.5 and 4000/3.0 (2), 3500/4.5 and 
3500/3.0 (3). Solid lines correspond to logg = 3.0, dashed ones to 
logg = 4.5, with a lithium abundance logiV(Li) = 3.5 



of the atmosphere. Note that the dependence of thermalization 
processes on the gravity is less pronounced. 

In the case of the subordinate lines we find that in the outer 
part of the atmosphere Si/B v > 1.0, while in the formation 
region of these lines (r ~ 0.01) is Si/B v ~ 1.0. The ratio 
Si/Bv for the A610.3 nm line is greater than for the A812.6 
nm line, because in the frequencies of the former line dB/dr 
is greater. 

3. 3. The profile of the resonance doublet 

Here we compare the LTE and non-LTE profiles of the reso- 
nance doublet. The absorption doublet is saturated in the case 
of stellar atmospheres with T cS < 4000 and logJV(Li) > 2.0. 
Thus, to describe the radiation transfer in the line core fre- 
quencies we extrapolated the Kurucz model atmospheres up 
to levels where line cores are formed, as described in Pavlenko 
(1984) and Magazzu et al. (1992). In Fig. | we show a few 
profiles computed for lithium abundance logiV(Li) = 3.5. We 
see that non-LTE cores are deeper in comparison with LTE, 
as found in Magazzu et al. (1992); see also Pavlenko (1992) 
and Carlsson et al. (1994). The Li I resonance line computed 
for model atmospheres with larger loggr shows stronger wings, 
both in LTE and non-LTE, due to the pronounced damping. 
However, we see very small differences in the non-LTE cores 
computed for model atmospheres with different logo's. In the 
case of unsaturated Li I resonance line the non-LTE profiles 
are shallower than the LTE ones (see Fig. 7 in Pavlenko et al. 
1995). 

The shape of the Li doublet becomes extremely important 
in some astrophysical contexts (e.g. Li isotopic ratio determi- 
nations). The difference between a LTE and a non-LTE line 
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Fig. 3. The LTE and non-LTE profiles of the Li I resonance dou- 
blet computed for a) 4000/3.0 (1) and 4000/4.5 (2) model atmo- 
spheres; b) 3500/3.0 (1) and 3500/4.5 (2) model atmospheres. Solid 
and dashed lines show LTE and non-LTE results, respectively. Dots 
indicate log g = 3.0 and open squares log g = 4.5. Note how the cores 
of non-LTE profiles do not differ for different logo's 



profile, calculated for the same equivalent width and for a 
model atmosphere 4000/3.0, is shown in Fig. ^. The largest 
differences are found in the core and core-wing transition part 
of the line. Note that these differences are negative in the core 
and positive in the wings. 

3.4- The curves of growth 
3.4.1. Resonance doublet 

The LTE and non-LTE curves of growth of the resonance 
doublet, computed for model atmospheres of several dwarfs 
and subgiants are presented in Table |l| In Fig. || we plot 
some of these curves, corresponding to the models 3500/3.0 
and 5000/4.5. We see that non-LTE abundance corrections 
A log N (Li) = log7V NLTE (Li) - log7V LTE (Li) are positive for 
unsaturated lines formed in the atmospheres of dwarfs and 
subgiants. As noted before, in that case the influence of ove- 
rionization of lithium dominates. The formation region of the 
non-LTE doublet shifts (with respect to LTE) towards deeper 
layers, where Si(T n \tc = 1) > B v (j\t c = 1), being r the optical 
depth at the frequency of the resonance line. A stronger dou- 
blet is formed where Si(r n i tc = 1) < B v (r\te = 1). As a result, 
the sign of AlogiV(Li) changes (Magazzu et al. 1992). Again, 
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Table 1. Curves of growth of the A670.8 nm line 
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Fig. 4. The difference between a LTE and a non-LTE Li I resonance 
doublet of the same equivalent width (W\ = 30.1 pm). Model at- 
mosphere is 4000/3.0 
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Fig. 5. The LTE and non-LTE curves of growth of the Li I reso- 
nance doublet computed for 3500/3.0 (1) and 5000/4.5 (2). Solid 
and dashed lines show the LTE and non-LTE results, respectively 



in the case of the subgiant 4000/3.0, when the doublet is ex- 
tremely strong AlogiV(Li) changes its sign once more. The 
corresponding line formation region moves towards the outer 
boundary of the atmosphere where the overionization increases 
with height and wins over the source function effect. This ten- 
dency can be seen also in the top panels of Figs. 16 and 17 in 
Carlsson et al. (1994) and in Fig. 4 in Pavlenko et al. (1995). 

Another result seems to be more interesting: our computa- 
tions show that non-LTE curves of growth computed for dwarfs 
and subgiants with Tig < 4000 K do not differ significantly in 
the case of unsaturated lines. In other words, whereas in LTE 
the equivalent width of the (unsaturated) Li I doublet is a func- 
tion of T e ff, loggr, and iV(Li), non-LTE computations give W\ 
as a function of only T e g and JV(Li). This is a very useful re- 
sult, because the determination of the gravity usually is more 
difficult than the assignment of an effective temperature. 




Fig. 6. The same as in Fig. [ij but for subordinate line at A610.3 nm 



3.4.2. Subordinate lines 

In Tables 2, 3 and Figs. 6, 7 we show the case of the subordinate 
Li I lines at A610.3 and A812.6nm. These lines are less intense 
than the resonance doublet and form in deeper parts of the 
atmosphere (Pavlenko et al. 1995). In spite of this, for these 
unsaturated lines we have also determined that the equivalent 
width in non-LTE practically does not depend on gravity for 
the model atmospheres with T c g < 4000K. 

The scarce dependence of all the considered Li I lines on 
log g may be explained in the following way. The statistical bal- 
ance of lithium in cool atmospheres is governed by the radiation 
field (Steenbock & Holweger 1984). The radiation field mean 
intensity crucially depends on the value of electron density (N c ) 
in the stellar photosphere. We found that differences of N e in 
the region 0.1 < n,, < 1.0 in the photosphere of stars with 
different logo's are much more pronounced for T e g > 4500 K 
than for lower T e g's. As a result, in the last case the intensity of 
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Table 2. Curves of growth of the A610.3 nm 


line 
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95.6 


87.3 


98.7 


89.0 


103.3 


94.7 




2.5 


2.3 


2.3 


2.5 


2.4 


2.6 


2.5 




3.0 


7.3 


7.1 


7.9 


7.4 


8.2 


7.7 


T cS = 6000A" 


3.5 


21.5 


20.9 


23.2 


22.0 


24.1 


22.8 




4.0 


56.8 


56.0 


60.5 


58.4 


62.9 


60.7 



the outgoing radiation field in the frequencies of lithium tran- 
sitions is less sensitive to logg. Moreover, the role of inelastic 
collisions in the formation of statistical equilibrium decreases 
with lowering T c s due to the decrease of N B . Note that the in- 
dependence from log g is seen only for unsaturated lines, which 
do not have extended, sensitive to logg, wings. 

3.4.3. Comparison with other works 

The LTE and non-LTE curves of growth presented here agree 
well with those used in the works of Martin et al. (1994) and 
Garcia Lopez et al. (1994), the differences in equivalent width 
being less than 5 percent. There is a good agreement also with 
Magazzu et al.'s (1992) results, except for cooler stars. Reasons 
and consequences of this disagreement at lower temperatures 
are discussed in Pavlenko et al. (1995). Note that in this pa- 
per we use a refined value of the damping constant for the Li I 



resonance doublet. We use also an extended opacity list, that 
makes the curves of growth presented here more reliable, es- 
pecially for lower effective temperatures. For these points see 
also Pavlenko et al. (1995). 

Carlsson et al. (1994) did not give the computed equiv- 
alent widths of lithium lines, but only non-LTE abundance 
corrections to apply to LTE curves of growth for the lithium 
resonance doublet. In order to compare our results with theirs, 
we use the data of Table 1 to compute our non-NLTE correc- 
tions. In Fig. ^| we present AS, i.e. the difference between our 
corrections (A us ) and those (Ac ar ) by Carlsson et al. (1994) 
for several models in common. 

We find that the agreement of these results is fairly 
good (see also Pavlenko 1995), the differences being less than 
0.1 dex. The greater differences are seen in strong Li I lines 
and are caused by the cumulative effect of different model at- 
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Table 3. Curves of growth of the A812.6 nm line 





logTV(Li) 






W x 


(mA) 










logs 


= 3.0 


log g 


= 3.5 


log g 


= 4.5 






LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 




0.5 


0.6 


0.4 


0.5 


0.4 


0.4 


0.4 




1.0 


1.8 


1.3 


1.5 


1.3 


1.2 


1.3 




1.5 


5.5 


4.1 


4.7 


4.1 


3.7 


4.2 


T of f = 3500/\ 


2.0 


16.8 


13.0 


14.5 


12.9 


11.5 


12.9 




2.5 


48.1 


38.7 


41.8 


38.3 


33.8 


38.2 




3.0 


117.3 


101.8 


104.7 


100.3 


87.2 


100.0 




3.5 


218.3 


210.8 


203.0 


208.3 


179.2 


207.4 




4.0 


314.3 


333.0 


302.4 


330.9 


284.5 


331.5 




1.0 


0.6 


0.4 


0.6 


0.4 


0.5 


0.4 




1.5 


2.1 


1.3 


2.0 


1.4 


1.6 


1.4 




2.0 


6.4 


4.2 


6.3 


4.3 


4.9 


4.4 


T eff = 40007\ 


2.5 


19.6 


13.0 


19.3 


13.4 


15.2 


13.6 




3.0 


55.3 


38.9 


54.6 


40.0 


43.9 


40.6 




3.5 


131.1 


103.2 


129.8 


105.8 


109.4 


107.3 




4.0 


232.1 


214.3 


230.5 


218.8 


209.8 


222.9 




2.0 


2.2 


1.4 


2.2 


1.5 


2.2 


1.6 




2.5 


6.8 


4.5 


6.8 


4.7 


6.9 


5.1 


T e ff = 4500if 


3.0 


20.6 


14.0 


20.7 


14.5 


20.9 


15.9 




3.5 


58.0 


41.5 


58.3 


43.0 


58.7 


47.1 




4.0 


136.3 


108.1 


136.7 


111.7 


137.9 


122.4 




2.0 


0.8 


0.6 


0.8 


0.6 


0.8 


0.7 




2.5 


2.5 


1.9 


2.6 


2.0 


2.7 


2.1 


T cff = 5000T\ 


3.0 


7.7 


5.9 


8.0 


6.1 


8.3 


6.7 




3.5 


23.2 


18.2 


24.0 


18.8 


25.0 


20.5 




4.0 


64.1 


52.6 


65.9 


54.4 


68.4 


59.4 




2.5 


1.0 


0.9 


1.1 


1.0 


1.2 


1.0 




3.0 


3.2 


3.0 


3.4 


3.1 


3.6 


3.3 


r eff = 55007\ 


3.5 


10.0 


9.3 


10.6 


9.6 


11.3 


10.2 




4.0 


29.7 


27.9 


31.4 


28.8 


33.1 


30.8 




2.5 


0.5 


0.5 


0.5 


0.6 


0.6 


0.6 




3.0 


1.5 


1.6 


1.7 


1.7 


1.8 


1.8 


T off = 60007^ 


3.5 


4.8 


5.1 


5.4 


5.5 


5.6 


5.6 




4.0 


14.8 


15.7 


16.5 


16.8 


17.0 


17.3 



mospheres, opacity sources, cross sections, model atoms. Note 
that the slope of curves of growth for strong lithium lines de- 
creases. We point out that the differences in the computed 
non-LTE abundance corrections are smaller than the correc- 
tions themselves. We also note that the data for weak lines 
formed deep in the atmosphere agree quite well. It is evident 
that these lines are formed in regions in the deep photosphere, 
where the temperature structures of models by different au- 
thors are practically the same. 

Our referee, Dr. Carlsson, kindly informed us about the re- 
sults of a comparison between our equivalent widths and those 
computed by Carlsson et al. (1994) for the 5000/4.5 model at- 
mosphere (solar metallicity). For LTE he found an agreement 
better than 10% with a difference of 5% for most abundances. 
The non-LTE equivalent widths agree within 15%, but typi- 
cally within 10%. 



4. Conclusions 

In this paper we give, in tabular form, the LTE & non-LTE 
curves of growth of lithium resonance (A670.8 nm) and subor- 
dinate lines (A610.3 and 812.6 nm) computed using the newest 
Kurucz (1993) models. These data may be used for lithium 
abundance determinations as well as for the redetermination of 
abundances obtained so far in the frame of the LTE approach. 
In particular, we point out the importance of the A812.6 nm 
data. This line lies in the red part of the spectrum and its 
blending is not as strong as for the other two lithium lines. 
The use of this line may give more accurate results. 

Our theoretical curves of growth cover a different effective 
temperature range than in Carlsson et al. (1994). In fact, our 
main interest is shifted towards lower temperatures which were 
not considered by these authors. We note a good agreement 
between our results and those obtained in other works, despite 
the differences in non-LTE procedure details. 



12 3 
log JV(Li) 



12 3 4 

log AT (Li) 



Fig. 7. The same as in Fig. H, but for subordinate line at A812.6 nm 



Let us note that lithium in the stellar atmospheres consid- 
ered in this paper exists mainly in ionized form. So the rel- 
atively small deviations from LTE obtained in our work for 
the coolest models of our grid have not to be considered as 
an obvious result. Indeed, with decreasing effective tempera- 
ture we have a dramatic drop of the electron pressure. On the 
other hand, opacities in the frequencies of lithium transitions 
increase due to the strong molecular absorption. In the case of 
high lithium abundances the resonance doublet becomes very 
strong with saturated core and extended wings are formed deep 
in the atmosphere, very close to LTE conditions. In the case 
of T c ff < 3500 K lithium exists in stellar atmospheres mainly 
in the form of neutral atoms, so non-LTE effects should not be 
large a priori. This was confirmed by direct computations by 
Pavlenko et al. (1995). 

For the coolest stars of our grid (3500 < T cti < 4000 K) we 
find that equivalent widths and profiles of unsaturated Li I lines 
computed in non-LTE show only a very weak dependence on 
the log g parameter. Instead, we note that the corresponding 
LTE curves of growth show a quite strong dependence on grav- 
ity. The use of non-LTE curves of growth gives a chance to min- 
imize the possible effects of errors in log g in lithium abundance 
determinations using non-saturated absorption lines. Moreover 
we note that, using LTE calculations to perform Li abundance 
determinations, the knowledge of gravity seems not crucial for 
T e ff > 4500 K, but it appears important for cooler stars. How- 
ever, when using non-LTE calculations, the reverse holds. 

In this paper we have discussed results for classical model 
atmospheres. The cool star spectra often show emission lines 
created in chromospheric-like features (CLF). However, as 
shown in Pavlenko (1995) and Pavlenko et al. (1995), CLF 
affect mainly the LTE results (both profiles and curves of 
growth). On the contrary, non-LTE results show weak depen- 
dence on CLF. Still in the case of the strongest CLF producing 
veiling the lithium lines may be severely affected. Let us note 



Fig. 8. The difference of non-LTE abundances corrections AS = A UH 
- Acar, computed by us and Carlsson et al. (1994), respectively, 
for dwarf model atmospheres 4500/4.5 (dots), 5000/4.5 (squares), 
5500/4.5 (crosses), 6000/4.5 (stars) 



that 

— only the youngest and the most active stars may have such 
CLF producing additional continuum in the visible part of 
the spectrum; 

— these results may be interpreted as a weak dependence of 
the (non-LTE) lithium lines on the structure of the outer- 
most layers of model atmospheres. At the present time we 
cannot be sure that we describe these layers properly in 
the frame of the classical one-dimensional approach. 

These two reasons give additional support to the use of 
non-LTE results for numerical analysis of lithium lines in stellar 
spectra. 

Acknowledgements. We thank Dr. M. Carlsson, our referee, for 
valuable comments which improved this paper. 
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